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DELTA-SIGMA  UHF  DIGITAL  WAVEFORM  GENERATOR 


1.  INTRODUCTION 

The  Naval  Research  Laboratory  (NRL)  has  developed  and  built  a  prototype  Digital  Waveform 
Generator  (DWG)  in  the  ultra-high-frequency  (UHF)  range.  The  objective  of  this  project  was  to  develop, 
test,  document,  and  deliver  a  preproduction-ready  high-performance  prototype  DWG.  This  prototype  uses 
Delta-Sigma  (A-S)  modulation  techniques,  which  permit  arbitrary  and  accurate  waveform  generation,  to 
meet  the  demanding  and  diverse  waveform  requirements  of  future  radar  applications,  including  linear-FM 
(LFM)  and  continuous  wave  (CW)  signals.  This  DWG  also  allows  for  the  generation  of  waveforms  at 
other  frequencies  by  upconversion  or  downconversion.  The  work  described  here  is  based  on  some  of  the 
earlier  work  of  Scholnik  [1]  on  A-S  topology  and  Harrell  [2]  on  single-bit  digital  signal  synthesizers. 
Some  of  the  efforts  on  this  project  were  previously  reported  in  Reference  3. 

The  development  of  this  DWG  concept  was  relatively  simple  and  required  only  minimal  computing 
effort.  It  provides  predicted  low  phase  and  spurious  noise.  Much  of  the  design  uses  field  programmable 
gate  arrays  (FPGAs)  for  single-bit  digital  waveform  generation.  Four  filter  topologies  were  initially 
considered,  including  cascade,  hybrid,  parallel,  and  transposed.  However,  the  cascade  and  parallel  forms 
were  eliminated  because  they  imposed  a  heavy  computational  burden  on  the  system.  Following  analysis 
of  the  transposed  topology,  quantization  error  in  higher-order  filters  lead  to  the  selection  of  the  hybrid 
form  because  it  provides  the  best  performance.  A  12*  order  hybrid  filter  was  selected  and  implemented 
using  FPGAs.  FPGAs  are  rapidly  improving  with  embedded  multipliers,  more  memory  availability  within 
the  chip  architecture,  and  greater  ease  of  use.  Utilizing  this  advancing  technology,  NRL  coded  a  FPGA  to 
generate  a  single-bit  A-S  waveform. 

The  NRL  development  demonstrates  that  a  simple-to-code  single-bit  A-S  DWG  can  provide  the  same 
resolution  as  a  16-bit  or  greater  digital-to-analog  converter  (DAC)  DWG.  Based  on  the  results  of  the 
effort  described  here,  this  technology  shows  the  promise  of  providing  the  diverse  waveform  and  improved 
accuracy  requirements  of  future  high-performance  Navy  radar  applications  with  reduced  complexity  at 
relatively  low  cost. 

2.  DELTA-SIGMA  ALGORITHM 

The  A-S  filter  algorithm  or  A-S  modulator  is  used  to  generate  a  single-bit  waveform  from  the  digital 
input  data  with  a  CW  center  frequency  between  406  and  450  MHz,  both  with  and  without  a  4-MHz 
bandwidth  LFM,  with  a  nominal  pulsewidth  of  10  ps  and  an  available  filter  bandwidth  of  50  MHz. 
Initially,  four  different  A-S  digital  filter  topologies  were  considered:  (1)  cascade;  (2)  hybrid;  (3)  parallel; 
and  (4)  transposed.  The  order  M  of  a  A-S  modulator  is  indicated  by  the  number  of  feedback  loops.  The 
cascade  and  parallel  forms  allow  for  a  growing  critical  path  in  proportion  to  increasing  filter  order.  They 
are  thus  computationally  more  intensive  and  challenging  for  the  FPGA  hardware  and  were  therefore 
eliminated  from  further  consideration.  Only  the  transposed  and  hybrid  forms  were  considered  here  in 
order  to  achieve  the  intended  >  90  dB  dynamic  range  for  CW  signals. 
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The  basic  operation  of  the  A-S  modulator  in  encoding  waveforms  can  be  understood  more  intuitively 
by  demonstration  [4].  A  simple,  first-order  modulator  is  shown  in  Fig.  1.  The  output  of  the  comparator, 
V(z),  is  also  the  output  of  the  loop,  which  is  full-scale  ±1  V.  At  the  summation,  either  +1  V  or  -1  V  is 
added  to  the  input  voltage.  The  result  is  the  input  to  the  A-S  digital  filter.  The  digital  filter  can  implement 
the  filter  equation  using  the  two  fundamental  arithmetic  operations,  multiplication  and  addition  (or 
accumulation).  If  the  output  of  the  digital  filter,  Y(z),  is  greater  than  0  V,  V(z)  becomes  +1  V;  if  it  is  less 
than  0  V,  V(z)  becomes  -1  V.  A  no-change  condition,  however,  causes  the  modulated  signal  to  remain  at 
the  same  ±1  V  state  of  the  previous  sample.  Each  operation  occurs  once  during  each  clock  cycle.  The 
output  of  the  loop,  either  +1  V  or  -1  V,  will  be  saved  in  a  single-bit  data  file  as  either  a  1  or  0, 
respectively.  Then  the  digital  single-bit  data  series  is  converted  to  an  analog  signal  by  using  a  high-speed 
single-bit  DAC.  In  general,  first-  and  second-order  modulators  are  stable,  but  stability  must  be  considered 
in  higher-order  units. 


Fig.  1  — A  typical  A-S  modulator  to  encode  digital  waveforms 


2.1  Delta-Sigma  Filter  Topologies 

Both  transposed  and  hybrid  A-S  approaches  were  considered  as  described  in  the  next  two  sections  of 
this  report.  The  approaches  included  a  comparison  of  transposed  filters  of  various  orders.  The  transposed 
filter  quantization  error  in  higher-order  (e.g.,  M  =  8)  filters  was  considered  and  hybrid  filters  of  various 
orders  were  compared. 

2.1.1  Transposed  Form 

The  transposed  filter  has  a  short  critical  path  of  two  full  additions  and  one  multiply.  It  thus  has  the 
advantage  of  a  short  critical  path  and  an  efficient  finite-impulse  response  (FIR)  implementation.  It  is, 
however,  extremely  sensitive  to  coefficient  quantization  accuracy.  Transposed  filters  were  analyzed  for 
three  orders  (M  =  2,  4,  8).  In  particular,  the  sensitivity  of  quantization  error  vs  signal-to-noise  ratio  (SNR) 
was  investigated  and  is  addressed  below. 

Transposed  filters  of  order  M  are  implemented  with  an  efficient  FIR  filter,  where  two  full  adders  and 
one  multiplier  are  used  in  each  stage.  For  an  N-length  discretization  [X(l),  X(2),  ...,  X(N)]^  of  the  input 
signal  X(z),  the  N  equations  for  the  M*  order  transposed  A-S  filter  are  given  by 

M 

Y(n)  ==  ^  {amY (n  -  m)  +  bm[X(n  -  m)  +  V(n  -  m)] } ,  for  n  >  m, 

m=l 


where  am  and  bm  are  the  m^^  coefficients  of  the  filter,  X(n)  is  the  input  signal  at  the  n^^  discrete  value  of  z 
for  n  =  1,. .  .,N,  and  V(n)  is  the  output  of  the  A-S  filter  at  the  n^^  discrete  value  of  z.  In  particular. 
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V(n)  = 


J+  l,if  Y(n)>0 
Y(n)<0. 


A  typical  transposed  A-S  digital  filter  for  M  =  4  is  shown  in  Fig.  2.  The  associated  transposed-filter 
coefficients  for  M  =  2,  4,  8,  as  shown  in  Table  1,  were  computed  using  the  MatLab  A-S  coefficient- 
generator  subroutine.  The  A-S  digital  filter  provides  noise  shaping  within  the  filter  bandwidth.  The  CW 
spectra  of  the  transposed  filter  outputs  are  shown  in  Fig.  3  for  (a)  M  =  2,  (b)  M  =  4,  and  (c)  M  =  8.  On 
observing  Fig.  3,  it  can  be  seen  that  as  the  filter  order  increases  the  SNR  increases  within  the  filter 
bandwidth,  which  is  50  MHz  in  this  example.  However,  with  higher-order  filters  such  as  M  =  8, 
coefficient  quantization  decreases  the  SNR  within  the  filter  bandwidth. 


Fig.  2  —  A  typical  transposed  A-S  digital  filter  for  M  =  4 


Table  1  —  Transposed  Filter  Coefficients  for  M  =  2,  4,  8 


Filter  Order  M 

Coefficients  ai,  aM 

Coefficients  bi,  . . bM 

9 

ai  =  0.3956 

bi  =  0.9503 

Z 

a2  =  0.5878 

b2  =  -l 

ai  =  0.3988 

bi  =  1.8988 

A 

a2  = -1.0661 

b2  =  -2.8973 

4 

a3  =  0.8812 

b3  =  1.8988 

a4  =  -0.5567 

b4  =  -l 

ai  =0.3851 

bi  =  3.7978 

a2  = -1.8220 

b2  =  -9.3983 

a3  =  4.0464 

b3  =  14.7971 

8 

a4  =  -6.0648 

b4  = -17.5997 

as  =  6.0185 

bs  =  14.7971 

ae  =  -4.3489 

be  =-9.3983 

ay  =  1.9440 

by  =  3.7978 

a8= -0.5576 

b8=  -1 
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(C) 


Fig.  3  —  CW  spectra  using  A-S  transposed  filters:  (a)  M  =  2,  (b)  M  =  4,  and  (c)  M  =  8. 
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To  show  that  SNR  reduction  in  a  higher-order  filter  is  sensitive  to  coefficient  quantization,  hi  = 
3.7978  for  the  filter  spectra  of  Fig.  3(c)  for  M  =  8,  as  shown  again  for  comparison  in  Fig.  4(a),  is  changed 
to  hi  =  3.7976  with  the  resulting  spectra  shown  in  Fig.  4(b).  On  comparing  Figs.  4(a)  and  4(b),  it  can  be 
observed  that  this  filter  is  extremely  sensitive  to  coefficient  quantization  for  higher-order  filters.  In  this 
example,  the  sidelobe  levels  increase  by  at  least  20  dB.  The  A-S  hybrid  filter,  which  has  the  same  critical 
path  as  the  transposed  filter,  is  considered  in  the  following  section. 


(a) 


(b) 


Fig.  4  —  CW  spectra  using  a  A-Z  transposed  filter  for  M  =  8  with  (a)  bi  =  3.7978  and  (b)  bi  =  3.7976 
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2.1.2  Hybrid  Form 

The  hybrid  filter  has  the  same  short  critical  path  as  the  transposed  filter  with  two  full  additions  and 
one  multiply.  It  thus  also  has  the  advantage  of  a  short  critical  path,  but  poles  on  the  unit  circle  reduce  the 
multiply  count.  However,  the  zeros  may  be  sensitive  to  coefficient  quantization  accuracy.  The  impact  of 
hybrid  filters  is  determined  and  analyzed  for  M  =  2,  4,  8,  12,  16.  The  equations  for  the  order  hybrid 
A-S  filter  are  given  by 

YM(n)  =  aM  [YM.i(n-l)]  +  bM  [X(n-l)  +  V(n-l)]; 

YM-k  (n)  =  aM-k[YM-k(n-l)]  +  bM-k  [X(n-l)  +  V(n-l)]  +  YM.k+i(n-l);  for  k  =  1,3,...,M-1 

YM-k  (n)  =  aM-k[YM-k-i(n-l)]  +  bM-k  [X(n-l)  +  V(n-l)]  +  YM-k+i(n-l);  for  k  =  2,4,...,M-2 


where  aM  and  bM  are  the  coefficients  of  the  filter,  X(n)  is  the  input  signal  at  the  n^^  discrete  value  of  z 
for  n  =  1,. .  .,N,  and  V(n)  is  the  output  of  the  A-S  filter  at  the  n^^  discrete  value  of  z.  In  particular, 

r+I,  ifY(n)>0 

V(n)= 

[-I,ifY(n)<0. 

A  typical  hybrid  A-S  digital  filter  for  M  =  4  is  shown  in  Fig.  5.  The  associated  hybrid  filter 
coefficients  for  M  =  2,  8,  12,  16,  as  shown  in  Table  2,  were  computed  using  the  MatLab  A-S  coefficient- 
generator  subroutine.  The  CW  spectra  of  the  hybrid  A-S  filter  outputs  are  shown  in  Fig.  6  for  (a)  M  =  2, 
(b)  M  =  8,  (c)  M  =  12,  and  (d)  M  =  16  A-S  hybrid  filters.  The  corresponding  equations  for  the  4*  order 
hybrid  filter  can  be  easily  written  as 

Y4(n)  =  a4  [Y3(n-1)]  +  b4  [X(n-l)  +  V(n-l)]; 

Y3(n)  =  a3  [Y3(n-1)]  +  b3  [X(n-l)  +  V(n-l)]  +  Y4(n-1); 

Y2(n)  =  a2  [Yi(n-l)]  +  b2  [X(n-l)  +  V(n-l)]  +  Y3(n-1); 

Yi(n)  =  ai  [Yi(n-l)]  +  bi  [X(n-l)  +  V(n-l)]  +  Y2(n-1). 


Fig.  5  —  A  typical  hybrid  form  A-S  digital  filter  for  M  =  4 
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Table  2  —  Hybrid  Filter  Coefficients  for  M  =  2,  8,  12,  16 


Filter  Order  (M) 


Coefficients  ai,....,  Em 


Coefficients  b 


ai  =  0.9503 

bi  =  0.3956 

a2  =  -l 

by  =  -0.5878 

ai  =  1.0440 

bi  =  0.3856 

a2  =  -1 

by  =  -0.7613 

a3  =  0.9876 

by  =  -0.1778 

a4  =  -1 

b4  =  0.2445 

as  =  0.9125 

bs  =  0.0375 

a6  =  -l 

be  =  -0.0442 

ay  =  0.8537 

by  =  -0.0035 

a8  =  -l 

bs  =  0.0034 

ai  =  1.0517 

bi  =  0.3856 

a2  =  -l 

by  =  -0.7686 

a3  =  1.0225 

by  =  -0.1749 

a4  =  -1 

b4  =  0.2583 

as  =  0.9765 

bs  =  0.0389 

a6  =  -l 

be  =  -0.0588 

ay  =  0.9238 

by  =  -0.0051 

a8  =  -l 

bs  =  0.0090 

as  =  0.8764 

bs  =  0.0004 

aio  =  -l 

bio  =  -0.0008 

ai  \  —  0.8456 

bii  =0.0000 

ai2  =  -1 

bi2  =  0.0000 

ai  =  1.0547 

bi  =  0.3369 

a2  =  -l 

b2  = -0.6781 

as  =  1.0371 

b3  =  -0.1191 

a4  =  -1 

b4  =  0.1997 

as  =  1.0078 

b5  =  0.0182 

a6  =  -l 

be  =  -0.0421 

ay  =  0.9705 

by  =  -0.0006 

a8  =  -l 

bs  =  0.0062 

as  =  0.9299 

bs  =  -0.0003 

aio  =  -l 

bio  =  -0.0005 

an  =  0.8916 

bii  =0.0000 

ai2  =  -1 

bi2  =  0.0000 

ai3  =  0.8610 

bi3  =  0.0000 

ai4  =  -l 

bi4  =  0.0000 

ais  =  0.8424 
ai6  =  -l 

bis  =  0.0000 
bie  =  0.0000 
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Fig.  6  —  CW  spectra  using  A-Z  hybrid  filters:  (a)  M  =  2,  (b)  M  =  8,  (c)  M  =  12,  and  (d)  M  =  16 

(see  next  page  for  (c)  and  (d)) 
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Frequency  (MHz) 
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Fig.  6  —  (continued)  CW  spectra  using  A-Z  hybrid  filters:  (a)  M  =  2,  (b)  M  =  8,  (c)  M  =  12,  and 

(d)M=  16 


Figure  6  indicates  that  the  dynamic  range  of  the  signal  increases  with  filter  order.  However,  with  the 
use  of  higher-order  filters  such  as  with  M  =  16,  coefficient  quantization  reduces  the  SNR  to  an 
unacceptable  level.  It  should  be  noted  that  the  93  dB  SNR  shown  in  Fig.  6(c)  is  comparable  to  a  16-bit  or 
greater  conventional  DAC  DWG.  These  figures  also  show  the  spectrum  of  a  CW  signal  at  425  MHz  with 
more  than  90  dB  desirable  dynamic  range  using  a  A-S  hybrid  filter  for  M  =  12.  Based  on  these  simulated 
filter  outputs,  the  hybrid  digital  filter  with  M  =  12  is  chosen  to  generate  UHF  waveforms  in  the  range  of 
406  to  450  MHz. 
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3.  SYSTEM  HARDWARE  DESCRIPTION 

The  DWG  is  designed  to  take  advantage  of  current  technology  advances  in  digital  circuitry  including 
the  use  of  FPGAs.  An  FPGA  system  provides  a  flexible  programmable  computer-logic  environment  in 
which  to  generate  waveforms  using  A-S  modulation  techniques.  The  general  system  diagram  is  shown  in 
Fig.  7.  This  section  describes  the  high-level  functionality  of  the  UHF  DWG  in  implementing  A-S 
techniques. 

A  MatLab  code  was  used  to  generate  the  digital  waveform  at  the  UHF  frequency  range  of  406  to  450 
MHz  with  an  oversampling  frequency  of  2488  MHz.  Typical  waveform  parameters  include:  (1)  a  UHF 
frequency  of  425  MHz;  (2)  a  bandwidth  of  4  MHz;  (3)  a  pulsewidth  of  10  ps;  and  (4)  a  pulse  repetition 
interval  (PRI)  of  0.1  to  10  kHz.  This  digital  waveform  code  is  loaded  onto  BlockRAM  memory  contained 
in  the  FPGA  chip.  The  waveform  is  then  filtered  using  a  A-S  algorithm.  The  output  of  this  filter  is  a 
single-bit  A-S  waveform.  The  serial  output  of  the  single-bit  waveform  is  converted  from  a  single-bit  to  a 
20-bit  parallel  format  and  is  loaded  in  a  selected  BlockRAM  memory  location  to  be  read  out  to  the 
synchronizer.  The  functionality  of  these  FPGAs  is  shown  in  Fig.  8. 


Fig.  7  — DWG  board 


Fig.  8  —  FPGA  functionality 


4.  SIMULATION 

Initial  simulation  of  the  A-S  algorithm  is  conducted  in  a  MatLab  environment.  Following  successive 
MatLab  runs,  the  desired  filter  characteristics  are  programmed  into  the  FPGA  environment. 

4.1  MatLab  Simulation 

The  procedure  for  accomplishing  the  MatLab  simulation  is:  (1)  generate  the  coefficients  for  the 
chosen  filter  using  the  MatLab  A-S  coefficient-generator  subroutine;  (2)  write  the  mathematical  code  for 
the  selected  A-S  digital  filter;  (3)  generate  digital  data  representative  of  a  signal  with  a  center  frequency 
from  406  to  450  MHz  supporting  a  4  MHz  bandwidth  LFM  modulation  format  with  a  nominal  pulsewidth 
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of  10  |Lis;  and  (4)  output  a  single-bit  A-S  waveform  from  the  A-S  filter.  For  a  A-S  hybrid  filter  of  M  =  12, 
the  time-domain  modular  outputs  and  the  resulting  spectra  from  a  CW  signal  and  a  4-MHz  bandwidth 
LFM,  both  with  nominal  pulsewidth  of  10  ps,  are  shown  in  Figs.  9  and  10.  The  spectra  are  generated  as 
20  logio  of  the  fast  Fourier  transform  (FFT)  of  the  time-domain  outputs.  After  analyzing  these  figures,  one 
can  conclude  that  the  CW  signal  at  the  filter  frequency  range  of  406  to  450  MHz  could  produce  an  output 
with  more  than  a  90  dB  dynamic  range.  In  the  case  of  LFM,  about  a  46  dB  dynamic  range  is  expected 
because  of  the  accumulation  of  energy  loss  due  to  the  4-MHz  bandwidth  of  the  signal.  As  expected,  the 
spectrum  is  not  well  constrained  in  frequency  beyond  the  50-MHz  band-limited  bandwidth. 

4.2  FPGA  Simulation 

A  similar  procedure  provides  a  simulation  in  the  FPGA  environment.  The  spectra  resulting  from  a 
CW  signal  and  a  4  MHz  bandwidth  LFM  using  a  A-S  hybrid  filter  for  M  =  12,  both  with  nominal 
pulsewidth  of  10  ps,  are  shown  in  Fig.  11.  On  comparing  Figs.  9,  10,  and  11,  it  can  be  concluded  that 
similar  filter  output  dynamic  ranges  are  attainable  in  the  FPGA  environment. 

5.  HARDWARE  IMPLEMENTATION 

Data  generated  using  the  A-S  algorithm  were  loaded  into  the  NRL  UHF  A-S  waveform  generator  via 
an  Agilent  pattern  generator  as  well  as  by  means  of  a  Virtex-II  Pro-P7  FPGA  development  board. 
Waveforms  produced  by  the  hardware  implementation  were  nearly  as  high  in  quality  as  the  spectra  of  the 
waveforms  shown  in  Fig.  11.  In  the  FPGA  environment,  the  Xilinx  Rocket  I/O  architecture  technology 
was  chosen  to  sample  a  single-bit  A-S  output  at  a  data  rate  of  2.5  Gbps. 

6.  SUMMARY  AND  CONCLUSIONS 

The  Naval  Research  Laboratory  has  developed  and  built  a  prototype  UHF  Digital  Waveform 
Generator  based  on  the  Delta-Sigma  algorithm,  which  allows  arbitrary  waveform  generation.  The 
development  of  this  DWG  concept  was  relatively  simple  and  required  only  minimal  computing  effort.  It 
provides  predicted  low  phase  and  spurious  noise.  Much  of  the  design  uses  field  programmable  gate  arrays 
for  single-bit  digital  waveform  generation.  Four  filter  topologies  were  initially  considered,  including 
cascade,  hybrid,  parallel,  and  transposed.  However,  the  cascade  and  parallel  forms  were  eliminated 
because  they  imposed  a  heavy  computational  burden  on  the  system.  After  analyzing  the  transposed 
topology,  quantization  error  in  higher-order  filters  led  to  the  selection  of  the  hybrid  form  of  the  digital 
filter  because  it  performed  well.  A  12^^-order  hybrid  filter  was  selected  and  implemented  using  FPGAs. 
The  NRL  development  demonstrates  that  a  simple-to-code  single-bit  A-S  DWG  can  provide  the  same 
resolution  as  a  16-bit  or  greater  digital-to-analog  converter  DWG.  Based  on  the  results  of  the  efforts 
described  here,  this  technology  shows  the  promise  of  providing  the  diverse  waveform  requirements  of 
future  high-performance  Navy  radar  applications  at  relatively  low  cost. 
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(a) 


(b) 


Fig.  9  —  CW  time-domain  output  (a)  and  spectrum  (b)  at  425  MHz  using  a  MatLab  A-S  hybrid  filter  for  M  =  12 
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(a) 


Fig.  10  —  LFM  time-domain  output  (a)  and  spectrum  (b)  at  425  MHz  using  a  MatLab  A-Z  hybrid  filter  for  M  =  12 
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(a) 


Fig.  11  —  Spectra  of  (a)  CW  and  (b)  LFM  signals  at  425  MHz  using  an  FPGA  A-S  hybrid  filter  for  M  =  12 
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